Measurement of Differences Between J/ip and ip' Suppression in p-A Collisions 
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Measurements of the suppression of the yield per nucleon of J/rp and ip' production for 800 GeV/c 
protons incident on heavy relative to light nuclear targets have been made with very broad coverage 
in xf and pr- The observed suppression is smallest at xf values of 0.25 and below and increases at 
larger values of xf- It is also strongest at small pr- Substantial differences between the ip' and J/ip 
are observed for the first time in p-A collisions. The suppression for the if)' is stronger than that for 
the J/ip for xf near zero, but becomes comparable to that for the J/ip for ip > 0.6. 
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Strong suppression of the yield per nucleon of heavy 
vector mesons produced in heavy nuclei relative to that 
in light nuclei has been observed in proton and pion- 
nucleus collisions Similar effects have also been 

observed in heavy-ion collisions 0. This suppression 
exhibits strong kinematic dependences, especially with 
Feynman-x (xf) and transverse momentum (jot) of the 
produced vector meson. Since the suppression of heavy 
vector meson production in heavy-ion collisions is pre- 
dicted to be an important signature for the formation of 
the quark-gluon plasma (QGP), it is important to un- 
derstand the mechanisms that can produce similar ef- 
fects in the absence of a QGP. These mechanisms can be 
studied in proton-nucleus production of vector mesons 
where no QGP is presumed to occur. Many effects have 
been considered |p|-|lT| in attempting to describe the ob- 
served proton-induced charmonium yields from nuclear 
targets, e.g. absorption, parton energy loss, shadow- 
ing and feed-down from higher mass resonances, but it 
is clear that no adequate understanding of the problem 
has been achieved. Even the absolute cross sections are 
poorly understood due to poor knowledge of the produc- 
tion mechanism, and most models ignore or use naive pic- 
tures of the space-time evolution of the cc pair. Recogniz- 
ing that the production and suppression mechanisms can 
be identified by their strong kinematic dependences, it is 



crucial to have new data with broad kinematic coverage 
to challenge comprehensive descriptions of charmonium 
production in nuclei. 

Here we report new high statistics measurements made 
in Fermilab E866/NuSea of the nuclear dependence of 
J/ip and ip' production for proton- nucleus collisions on 
Be, Fe, and W targets. Over 3 x 10 6 J/ip's and 10 5 ip n s 
with xf between —0.10 and 0.93 and pt up to 4 GeV/c 
were observed. Previous measurements in E772 0j and 
E789 have suffered from limited pt acceptance and 
limited statistics at larger values of xf , both of which 
are greatly extended in these new data. 

E866/NuSea used a 3-dipole magnet pair spectrome- 
ter employed in previous experiments (E605 |12], E772, 
and E789), modified by the addition of new drift cham- 
bers and hodoscopes with larger acceptance at the first 
tracking station and a new trigger system (l3| . This 
spectrometer was also used for other measurements in 
E866/NuSea ||@. An 800 GeV/c extracted proton 
beam of up to 6 x 10 11 protons per 20 s spill bombarded 
the targets used in these measurements. A rotating wheel 
which was located upstream of either the first or second 
magnet held thin solid targets of Be, Fe and W with 
thicknesses corresponding to between 3% and 19% of an 
interaction length. After passing through the target, the 
remaining beam was absorbed in a copper beam dump 
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FIG. 1. Fit to the mass spectrum for the, 3e, target 
in the x F range from 0.00 to 0.05. Components in the 
fit are the J/ip, the t[>', Drell-Yan (long-dashed), randoms 
(short-dashed), and open charm (dot-dash). The solid curve 
represents the total of all fitted lineshapes and the dotted 
curve represents the continuum which is the sum of the 
Drell-Yan, randoms and open charm. 

located inside the large second magnet. Following the 
beam dump was a 13.4 interaction length absorber wall 
which filled the entire aperture of the magnet, eliminated 
hadrons, and assured that only muons traversed the spec- 
trometer's detectors. These muons were then tracked 
through a series of detector stations composed of drift 
chambers, hodoscopes and proportional tubes. Because 
of improvements in the trigger system, the coverage in pp 
was much broader than in previous experiments with this 
spectrometer (e.g. E772), extending to over 4 GeV/c. 
Beam intensity was monitored using secondary-emission 
detectors. 

Three magnetic field and target location configura- 
tions were used to span the full range in x F : small- 
xp (SXF, -0.1 < x F < 0.3), intermediate-XF (IXF, 
0.2 < x F < 0.6) and large-x^ (LXF, 0.3 < x F < 0.93). 
Detailed Monte Carlo simulations of the J/ip and ip' 
peaks and of the Drell-Yan continuum were used to gen- 
erate lineshapes in each bin in xp or in pp. For the Drell- 
Yan calculations we use MRST H NLO with EKS98 @ 
shadowing corrections. The contribution to the contin- 
uum from semi-leptonic decay to muons of open charm 
pairs was estimated using PYTHIA [[l8| and a small cor- 
rection, less than half the statistical uncertainties, was 
made for it in the SXF data set; but for the larger xp data 
sets it is negligible and no corrections were made. In ad- 
dition, a detailed construction of random muon pairs us- 
ing single-muon events (which also provided a good fit to 
the like-sign muon mass spectra) was used to account for 
the smooth random background underneath the peaks. 
A maximum-likelihood method was used for fitting that 
took into account the statistical uncertainty of the data 



and of the Monte Carlo and randoms JljJ . Figure 1 shows 
a typical fit to a mass spectrum using these components. 
Since the rates in the various detectors were nearly equal 
for the different targets, a correction for rate-dependent 
inefficiencies was not necessary. 

We present our results in terms of a, where a is ob- 
tained by assuming the cross section dependence on nu- 
clear mass, A, to be of the form a a — on X A a , where <jn 
is the cross section on a nucleon. For the SXF data, a was 
obtained using Be and two different thickness W targets, 
while for the IXF and LXF data, Be, Fe and W targets 
were used. The SXF data from the two W targets ver- 
ified that no corrections for secondary production were 
necessary. The pp dependence of a is shown in Fig. 2, 
where we see essentially the same increase in a for all x p 
ranges for both the J /if) and the ip', as well as for the 
200 GeV/c NA3 data |]. This increase is characteristic 
of multiple scattering of the incident parton and of the 
nascent cc in the final state. Note that for the IXF data 
the pt acceptance is truncated at about 2 GeV/c because 
a more restrictive trigger was used. 
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FIG. 2. a versus pr for J/ip (solid circles) and ip' (open 
boxes) production by 800 GeV/c protons. Results are shown 
for the three data sets - SXF, IXF and LXF (see text) - 
which have (if) = 0.055, 0.308 and 0.480, respectively. Only 
statistical uncertainties are shown. An additional systematic 
uncertainty of 0.5% is not included. Also shown are the NA3 
results at 200 GeV/c whose xf range can be seen in Fig. 4. 
The solid curves represent the parameterization discussed in 
the text. 

Previous experiments such as E772 have had a limited 
acceptance in pp which varied with xp. Since the value 
of a depends strongly on pp this can cause a distortion of 
the apparent shape of a versus xp. The improvements 
in the E866/NuSea trigger allowed a much broader pp 
acceptance than in these earlier measurements. However, 
for the lowest values of xp at each spectrometer setting 
our pp acceptance still becomes somewhat restricted. For 
the results presented here we have corrected the values of 
a(x F ) using a detailed simulation of our acceptance and 
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a differential cross section shape versus pp derived from 
our data. 
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FIG. 3. a for the J/ip versus xf for the three different data 
sets (top) and for the J/ip and ip' after the data sets are com- 
bined (bottom). Values are corrected for the pr acceptance, 
as discussed in the text. These corrections (Aq) have a maxi- 
mum value of 0.06 and are shown using the right-hand vertical 
scale in the top panel. The relative systematic uncertainty be- 
tween a for the J/ip and ip' is estimated to be 0.003, while 
the absolute systematic uncertainty is 0.01 in a; neither is in- 
cluded here. The solid curve represents the parameterization 
discussed in the text. 

The resulting dependence ofaon sj? is shown in Fig. 3 
and listed in Table I. The systematic uncertainty of 1% in 
the corrected a is dominated by the pp acceptance cor- 
rection, a for the J/ip is largest at values of xf of 0.25 
and below but strongly decreases at larger values of xf- 
For the ip' a is smaller than for the J/ip for xf < 0.2, 
remains relatively constant up to xp of 0.5 (becoming 
slightly larger than for the J/ip) and then falls to values 
consistent with those for the J/ip for xf > 0.6. The sig- 
nificance of the overall J/ip, ip' difference for xf < 0.2 
is about 4 sigma with respect to the statistical and rel- 
ative systematic uncertainties. This difference is con- 
sistent with less accurate results obtained by NA38 for 
p-A at 450 GeV/c ||, but is inconsistent with the quoted 
NA38 result that also included the p-p and p-d data from 
NA51. Although slightly larger a values for the ■0' than 
for the J/ip can be seen near xf = 0.55, we should point 
out that if instead we emphasize the velocity of the cc 
and plot a versus rapidity, then the agreement is quite 
good in this region. The reduced a at small xp is also ev- 
ident in Fig. 2 where a for the ip' falls consistently below 
that for the J/ip at low pp for the SXF data set. 

Our results for the J/ip a can be represented for con- 
venience by the simple parameterizations shown as solid 
lines in Figs. 2 and 3: a{x F ) = 0.960(1 - 0.0519a; F - 



TABLE I. a versus xf |20| for the J/ip and ip' . a is defined 
by a a = cr]v x A a and is equal to one if there is no suppression 
and the cross section scales simply as the number of nucleons. 
The average momentum fraction of the struck parton, X2 pf| , 
and the center-of-mass rapidity, y C m, are also shown. An 
additional systematic uncertainty of 1% is not included here. 



{Xf) {ycm)j/j, {x2)j/4, Otj/j, (Ucm)^' (x 2 )^' 



lip' 



— 0.065 
— 0.019 
0.027 
0.075 
0.124 
0.173 
0.223 
0.277 
0.332 
0.381 
0.431 
0.481 
0.531 
0.582 
0.632 
0.682 
0.732 
0.781 
0.828 
0.873 
0.913 



-0.390 
-0.115 
0.161 
0.433 
0.680 
0.896 
1.091 
1.288 
1.427 
1.551 
1.663 
1.764 
1.858 
1.945 
2.026 
2.098 
2.166 
2.228 
2.286 
2.338 
2.383 



0.1192 
0.0902 
0.0679 
0.0511 
0.0395 
0.0316 
0.0262 
0.0213 
0.0182 
0.0160 
0.0142 
0.0128 
0.0117 
0.0107 
0.00984 
0.00916 
0.00855 
0.00804 
0.00760 
0.00723 
0.00698 



0.962(7) 
0.953(3) 
0.955(2) 
0.955(2) 
0.952(3) 
0.955(4) 
0.951(6) 
0.917(11) 
0.916(6) 
0.888(7) 
0.875(6) 
0.852(5) 
0.831(5) 
0.811(5) 
0.789(6) 
0.772(5) 
0.772(7) 
0.739(10) 
0.760(17) 
0.733(32) 
0.611(71) 



-0.344 
-0.104 
0.132 
0.369 
0.588 
0.785 
0.974 
1.144 
1.281 
1.401 
1.512 
1.614 
1.705 
1.791 
1.869 
1.942 
2.009 
2.071 



0.1346 
0.1056 
0.0828 
0.0645 
0.0513 
0.0418 
0.0347 
0.0293 
0.0253 
0.0223 
0.0199 
0.0179 
0.0163 
0.0150 
0.0138 
0.0129 
0.0120 
0.0113 



0.929(24) 
0.918(14) 
0.931(11) 
0.932(11) 
0.931(12) 
0.913(18) 
0.940(22) 
0.923(36) 
0.910(18) 
0.884(16) 
0.885(15) 
0.874(16) 
0.881(16) 
0.845(20) 
0.751(25) 
0.790(36) 
0.718(49) 
0.727(69) 



0.338x|), and a(p T ) = A;(l + 0.0604p T + 0.0107p T 2 ), 
where A t = 0.870, 0.840, 0.782 and 0.881 for the SXF, 
IXF and LXF datasets and for the NA3 data, respec- 
tively. 

A comparison of our results with earlier data from 
E772 at 800 GeV/c § and also with NA3 at 200 GeV/c 
Ji| is shown in Fig. 4. It illustrates that the suppres- 
sion seen for J/ip production scales with xp but not 
with Pjfi^ above about 90 GeV/c, which corresponds to 
xp > 0.05 for our data and to if > 0.4 for NA3. Also of 
interest in these figures is a comparison of our results with 
those of E772. At the small-ccF end of the E772 data their 
published results drop significantly below our results. As 
was discussed above, the E772 data have severe narrow- 
ing of the pp acceptance for their smallest Xf bins, and 
a large correction that could easily bring the E772 points 
into agreement with our data is expected. Similar argu- 
ments hold for the E789 J/ip data (not shown) at small 
to negative xf [[|, where we estimate about an 8% cor- 
rection which would bring those results into agreement 
with ours. On the other hand, the large xf results from 
E789 Q appear to be high by more than their systematic 
uncertainty of 2.5%. 

The suppression of J/ip production near xp — is usu- 
ally thought to be caused by absorption, the dissociation 
of the cc pair by interactions with the nucleus or comovers 
H into separate quarks that eventually hadronize into D 
mesons. This model is supported by both the increased 
suppression of the ip' that we observe near xf = and 
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FIG. 4. a versus if and versus PjfJ? for the J/ip from 
E866/NuSea (800 GeV/c) (solid circles) compared to E772 
(open diamonds) and NA3 (200 GeV/c) (open squares) show- 
ing the scaling with xf (bottom) and lack of scaling with 
Pjf* (top). 

the absence of suppression of D meson production in the 
same kinematic region [ p2[ . At small xf, the velocity of 
the cc pair is low enough that it may hadronize within 
the nucleus, so the larger ip' would be absorbed more 
strongly p|pO[|. However, the observed constancy of a 
for both the J/ip and the ip' up to y cm ~ 1 complicates 
this interpretation since these models predict that faster 
cc pairs above xp ~ 0.1 would experience similar absorp- 
tion, whether they eventually hadronize outside the nu- 
cleus into a J/ip or a ip' . At larger values of xf, above 0.3, 
our data does show similar suppression for the J /ip and 
ip' . Furthermore, if absorption by the nuclear medium is 
the dominant suppression mechanism, the effect should 
scale with Pj'kF, but Fig. 4 shows that scaling breaks 
down in the middle of the region where we observe a to 
be constant. 

Shadowing of the small- x target gluon distributions is 
also thought to play a role in the observed suppression, 
but current estimates J|,[l7j predict at most a few percent 
drop in a(xp), even at the largest xf values. Also, as is 
seen for our data (but not shown) and was seen previously 
H, there is a lack of scaling with X2, which is related to 
that shown above for T> L jf® since X2 oc \/-p L A® . This 
appears to rule out large contributions from shadowing. 
Our studies Jl4| show that for Drell-Yan the dominant 



nuclear effect is shadowing of the anti-quark distributions 
and that the energy-loss of the incident quark is small. 
Although the incoming gluon's energy loss is expected 
to be larger by a color factor of 9/4 and the additional 
energy loss of the outgoing cc may be as large as that of 
a gluon, we still expect the overall contribution of energy 
loss to be small for resonance production. 

In conclusion, we have presented new data for the 
suppression of J/ip and ip' production in heavy versus 
light nuclei for 800 GeV/c proton- nucleus collisions. The 
kinematic coverage in xp (—0.10 to 0.93) and pp (0 - 
4 GeV/c) and statistical accuracy surpass that of previ- 
ous experiments. Corrections are made to the data to 
account for the narrowing pp acceptance at the smaller 
values of xf- The largest value of a (integrated over pp) 
of about 0.95 is seen at xf near 0.25 and below with 
strongly falling values for larger xf- The most striking 
new result is that the suppression for the ip' is stronger 
than that for the J/ip at xf near zero. 
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